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ABSTRACT 

New GMRT observations of the five- component pulsar B 185 7-26 provide detailed in- 
sight into its pulse-sequence modulation phenomena for the first time. The outer conal 
components exhibit a 7.4-rotation-period, longitude-stationary modulation. Several 
lines of evidence indicate a carousel circulation time P3 of about 147 stellar rotations, 
characteristic of a pattern with 20 beamlets. The pulsar nulls some 20% of the time, 
usually for only a single pulse, and these nulls show no discernible order or periodicity. 
Finally, the pulsar's polarization-angle traverse raises interesting issues: if most of its 
emission is comprised of a single polarization mode, the full traverse exceeds 180°; or 
if both polarization modes are present, then the leading and the trailing halves of the 
profiles exhibit two different modes. In either case the rotating vector model fails to 
fit the polarization- angle traverse of the core component. 

Key words: miscellaneous - methods:MHD — plasmas — data analysis - pulsars: 
general, individual (B1857-26) — radiation mechanism: nonthermal - polarization. 



I. INTRODUCTION 

Pulsar B1857-26, though discovered early (Vaughn & Large 
1970), has heretofore been studied almost entirely by aver- 
age methods. One of the small group of pulsars with five- 
component (M) profiles, its aggregate polarization has been 
investigated over a broad frequency band (Hamilton et al 
1977; McCulloch et al 1978; Manchester et al 1980; Morris 
et al 1980; van Ommen et al 1997; & Gould & Lyne 1998), 
and no effort to interpret the form of pulsar beams fails to 
mention it (Backer 1976; Lyne & Manchester 1988; Rankin 
1983a, 1986, 1990, 1993a,b; Mitra & Deshpande 1999). Lit- 
tle has been learned, however, regarding its pulse-sequence 
(hereafter PS) behaviour. Apart from one historical effort to 
determine its null fraction (Ritchings 1976), only recently 
did the fluctuation-spectral analyses of Weltevrede et al 
(2006, 2007; hereafter WES, WSE) identify its 7-rotation- 
period (hereafter Pi) modulation. PSR B 1857-26 is often 
compared with other prominent M pulsars, but it is not 
yet known whether it exhibits profile moding like B 1237+ 25 
{e.g., Srostlik & Rankin 2005) or core-component phenom- 
ena such as those seen in B0329+54 (Mitra et al 2007). 

We have carried out a sensitive new 325-MHz polari- 
metric observation of B 1857-26 using the Giant Metre- Wave 
Radio Telescope (hereafter GMRT) in Maharashtra, and a 
typical 200-pulse segment is given in Figure [1] Already we 
see evidence of interesting subpulse modulation. In the re- 
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Figure 1. A 200-pulse Stokes-/ pulse sequence from pulsar 
B1857-26 (central panel). The average profile and integrated in- 
tensity are given in the bottom and side panels, respectively. Note 
that the trailing component exhibits a regular 7-rotation-period, 
longitude-stationary modulation that is also perceptible in the 
leading component. The intensity scale in this and subsequent 
figures in arbitrary. 
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mainder of this paper, we describe the results of a series of 
PS analyses designed to address a number of these questions 
remaining about this star's individual pulse behaviour. To 
our surprise we found that we were able to identify and mea- 
sure a long period cycle which almost certainly can be inter- 
preted as the rotation interval of an emission-beam carousel 
in the manner of that found earlier for B0943+10 (Desh- 
pande & Rankin 1999, 2001; hereafter DR99, DROl). Such 
cycles are thought to be driven by ExB forces within a 
pulsar's polar flux tube along the lines of the Ruderman & 
Sutherland (1975; hereafter R&S) theory; however, all those 
so far determined or reliably estimateqj are longer or much 
longer than predicted by the theory. §11 describes the GMRT 
observations, and §111 our efforts to measure and model its 
polarisation-angle traverse, so as to better determine the 
emission geometry. In §IV we analyse the properties of the 
nulls, and in §V we outline the dynamics of its core compo- 
nent. §VI gives a discussion the fluctuation-spectral analy- 
ses, and §VII the evidence for a tertiary modulation cycle. 
§VIII then discusses the configuration of the pulsar's sub- 
beam carousel, and §IX gives a brief discussion and summary 
of the results. 



II. OBSERVATIONS 

Pulse-sequence polarisation observations of pulsar B 1857-26 
were acquired using the Giant Meterwave Radio Telescope 
(GMRT) north of Pune, India at 325 MHz on 2004 August 
27. The GMRT is a multi-element aperture-synthesis tele- 
scope consisting of 30 antennas which can be configured as 
a single dish. The polarimetry discussed here combined the 
array signals coherently in the upper of the two 16-MHz 
'sidebands' in a manner identical to that described for pul- 
sar B0329+54 in Mitra et al (2007). The pulsar back-end 
computed both the auto- and cross-polarized power levels, 
which were then recorded at a sampling interval of 0.512 
msec. A suitable calibration procedure as described in Mi- 
tra et al (2005) was applied to the recorded observations to 
recover the calibrated Stokes parameters I, Q,U and V. The 
duration of the PS is 19.4 minutes or 1945 pulses. In addi- 
tion, a total-power observation at 610 MHz with the same 
bandwidth, sampling time and duration was carried out on 
2007 February 4. 



III. EMISSION GEOMETRY 

A histogram of the polarisation-angle (hereafter PA) be- 
haviour of the 325- MHz observation is given in Figure [2] 
The upper panel shows the star's five components in to- 
tal intensity (Stokes /), the total Unear polarisation (Stokes 
^ NvQ^ + ^^]) the circular polarization (Stokes V 
[=LH-RH]); whereas the lower panel gives the PA density 
twice for clarity. Several aspects of this diagram deserve close 
discussion. 

First, the very flat PA curves under the wings of the 

1 B0834+06: Asgekar & Deshpande 2005; Rankin & Wright 
2007a; B0809+74: van Leeuwen et al 2003; B0834-26: Gupta et al 
2004; B1133+16: Herfindal & Rankin 2007; J1819+1305: Rankin 
& Wright 2007b. 
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Figure 2. PA histograms for the entire 325-MHz PS. The upper 
panels show the aggregate total power, total linear and circular 
polarisation (LH-RH), and the lower panel gives the polarization- 
angle (hereafter PA) density. The PA values are plotted twice for 
clarity, and the two curves indicate negative-going RVM fits to the 
PA traverse. One RVM fit (dotted curve) assumes that both sides 
of the profile represent the PPM and, clearly, it fails to describe 
the > 180° extent of this apparent PA rotation. A second RVM 
fit assumes that different OPMs predominate in the leading and 
trailing parts of the profile. The longitude origin is taken near the 
zero-crossing point of the antisymmetric Stokes V signature. 



profile appear compatible with a positive (equatorward) 
sense of the sightline impact angle f3 as our attempts to 
fit a rotating- vector model (hereafter RVM: Radhakrishnan 
& Cooke 1969; Komesaroff 1970) to the traverse bear out. 
Gould & Lyne's (1998) polarimetry confirms this aspect of 
the pulsar's PA traverse at meter wavelengths. At 1.4 GHz, 
however, the traverse under the wings of the profile is not at 
all flat and even shows different slopes on the leading and 
trailing sides of the profile (see also Johnston et al 2005). 
Clearly, such behaviour is difficult or impossible to reconcile 
with the RVM. 

Second, no sensible RVM fit0 could be obtained to the 
PA traverse under the assumption that most of the star's 
emission, throughout the profile, stems from a single orthog- 
onal polarisation mode (hereafter 0PM) — here the putative 
primary polarization mode (hereafter PPM). All the pub- 
lished observations seem to confirm that the PA traverse 
rotates negatively (clockwise) in the center of the profile, 
but if most of the emission reflects the PPM, then the full 
PA traverse significantly exceeds 180° both at meter and 
centimeter wavelengths. While the RVM PA excursion can 

^ The RVM fits done here use the method and convention de- 
scribed by Everett & Weisberg (2001). Here, the goodness of the 
fits is assessed by the reduced chi-square, which in the present case 
yielded significantly large values. Also a and (3 are correlated up 
to 98%. 
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exceed 180° for inner (poleward) sightline traverses, it can- 
not do so for outer ones, thus some different interpretation is 
needed. This difficulty is illustrated by the fit (dotted curve) 
in the bottom panel of Fig. (2] 

We also tried to fit B1857-26's PA trajectory making 
the somewhat unsavory assumption that different OPMs are 
dominant in the leading and trailing wings of the pulse pro- 
file (as we know of no excellent example of this configuration 
in other pulsars). This RVM fit is indicated by the thinner 
full curve in Fig. [21 and note that it is hardly excellent. It 
results in a large (some 3°) displacement between the PA in- 
flection point and the zero-crossing longitude of the circular 
polarization. 

One characteristic of the profile possibly supporting this 
interpretation is the deep minimum in the total linear polar- 
ization just following the core component, which is promi- 
nent in every published observation. This latter fit results 
in values for the magnetic latitude a and sightline-impact 
angle [3 of 155 and -1.8°, respectively (using the Everett & 
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Figure 3. Longitude-resolved fluctuation spectra for pulsar 
B1857-26 at 325 (top left and right) and 610 (bottom left) MHz. 
The full spectra are shown in the main panels, the integral spectra 
(for the outer components only) in the lower ones, and the side 
panels give the average profiles (along with the dashed bound- 
aries for computing the integral spectra). Most intervals of the 
PSs exhibit a bright narrow feature (top left, full 1900 PS) at 
0.14 c/Pi or a P3 value of some 7.4 Pi, whereas others show 
disorderly, weaker responses (right, pulses 1-600) which include 
the 0.14 c/Pi feature. The first 600 pulses of the 610-MHz PS 
(bottom left) exhibit a strong low frequency periodicity as well 
as 'sidebands' adjacent to the primary modulation feature. Here 
overlapping ffts of length 256 were used. 



Weisberg [2001] convention) — or equivalently an outer tra- 
verse with a, /3 some 25, +1.8° 

Finally, B1857-26's absolute fiducial PA orientation was 
determined by Johnston et al (2005) and also discussed by 
Rankin (2007), so it is important to understand the geome- 
try of the pulsar's 0PM emission. The two RVM fits (as well 
as the further conjecture) above make very different assump- 
tions about this geometry, and they result in different offsets 
between the centres of the linear and circular polarisation. 
Further study at multiple frequencies is needed to resolve 
these issues; but, fortunately, we need not settle them fully 
for our present purpose. 

^ Even had these fits been fully successful, they would not close 
all questions about the pulsar's PA behaviour. We also see a weak 
track under the core component suggesting a positive PA traverse 
(where the average PA tends to follow on the core's trailing edge) 
with a slope of perhaps +40-50°/°. Recent studies {e.g.^ Srostlik 
& Rankin 2005; Mitra et al 2007) show that core linear polarisa- 
tion does not always follow the RVM, and here it is concentrated 
in two "spots" near the +/- circular maxima, rather than indi- 
cating any clear RVM track. Were this the correct geometrical 
interpretation, a, (3 would be some 27, +0.6°. 
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Figure 4. Modulation amplitude (bottom panel dashed curve) 
and phase (top panel) for the pulse 1101-1612 sub-PS in Fig. [3] 
(top left), shown as a function of the average-profile power (solid 
curve bottom panel). The central panel shows the contour levels 
for the Irf (zoomed in between 0.08 to 0.2 c/Pi) and the left 
panel shows the longitude-averaged fiuctuation spectrum. Note 
that virtually all of the fiuctuating power here falls under the 
outer conal component pair (I &; V) — and primarily V. 



IV. NULLING BEHAVIOUR 

B 1857-26 exhibits some 20di3% null pulses, substantially 
more than the 10zb2.5% indicated by Ritchings' (1976) very 
old analysis. We computed null histograms for the observa- 
tions similar to Redman et a/'s fig. 12. In both, the pulsar 
maintains a relatively constant brightness (so scintillation 
effects are minimal) and the noise level is such that about 
0.3 </> provides a plausible threshold for null identifica- 
tion. The intensity distributions of the pulsar signal and 
noise overlap so that weak pulses cannot be reliably dis- 
tinguished from nulls. Nonetheless, we see that the pulsar's 
nulls are very short, with a strong propensity for 1 P\ dura- 
tion, and a maximum length in our observations of 4 Pi. Its 
bursts are also short, many having only a single period, and 
a mean burst length of only 5 Pi . We find little evidence that 
the nulls are non-random or even weakly periodic (using the 
methods of Herfindal & Rankin 2007). However, the aver- 
age profile of the putative nulls (not shown) does indicate 
significant power at the longitude of the central component, 
whereas the partial profiles of pulses just before and just 
after nulls show no significant difference in form from the 
total average profile. 



V. CORE DYNAMICS 

The PS properties of core-component emission have not been 
well studied, and B 1857-26 provides a further important op- 



portunity to examine its characteristics. This pulsar's core 
tends to be quite regular in intensity, as is suggested by 
Fig. [1] and further indicated by its rather small modulation 
indices at both 92 and 21 cms (WES, WSE). The core fea- 
ture is not at all symmetrical (or of Gaussian form) with its 
slow rise and steeper fall, and its constituent subpulses ex- 
hibit neither the intensity-dependent form nor the longitude 
shift seen in pulsar B0329+54 {e.g., Mitra ei al 2007). Even 
its antisymmetric circular polarisation is rather consistent, 
showing very nearly the same phase in every pulse. Inter- 
estingly, the strongest single pulses seem always to fall near 
the core peak and are highly left-hand circularly polarised; 
whereas the leading right-hand circular is somewhat weaker 
and steadier (see also WSE's asymmetric variance profile). 
Also, we see a slight change in the total linear polarisation 
from a double form at low intensities to an emphasis on the 
trailing lobe in the strongest pulses. Again, the circularly po- 
larised core signature in B 1857-26 bears careful examination 
as it fails to be fully antisymmetric in just the same manner 
as the Stokes / profile fails to be symmetric. Nonetheless, 
its zero- crossing point falls very close to the center of the 
conal profile. 



VI. FLUCTUATION-SPECTRAL ANALYSES 

Figure [3] gives longitude-resolved fluctuation (hereafter Irf) 
spectra for three different segments of the two PSs. As ex- 
pected from Fig. [1] periodic modulation is seen in the outer 
conal components (I and V, and its contribution is shown 
in the integrated fluctuation spectra in the bottom pan- 
els) which often exhibits a strong feature near 0.14 cycle/Pi 
(hereafter c/Pi). The 325-MHz Irf (top left) shows that this 
narrow feature dominates throughout the entire PS. Atop 
the broader response, its power falls in a single bin of a 512- 
length FFT (not shown) with a frequency of 0.135±0.001 
c/Pi or 7.41di0.06 Pi. The corresponding harmonic-resolved 
fluctuation spectrum (hereafter hrf, see DROl for details; 
also not shown) indicates that the respective positive and 
negative responses represent mostly amplitude modulation 
as expected. Figure |4] shows the modulation amplitude and 
phase for the interval 1101-1612 PS. Virtually all the power 
with this periodicity falls under the outer conal components 
(not the inner pair), and the flat phase gradient in this re- 
gion indicates a longitude- stationary modulation — ^just as 
expected for a highly central sight line traverse. 

Other sections of the 325-MHz PS, however, show less 
regularity. The right-hand Irf plot of Fig. [3] gives a repre- 
sentative example. Here, the fluctuations produce several 
broad features around the primary one. We emphasize that 
the 7- Pi modulation is associated with the outer conal com- 
ponents and primarily the trailing one; it cannot even be 
detected in the inner conal region of the proflle. The fluc- 
tuation spectra of the outer components is conflated in the 
integral spectra, but separate analyses (not shown) of the 
conal and core parts of the proflle show that the core con- 
tributes little (even noise power) to the integral spectra. In 
summary, only the outer cone fluctuates signiflcantly at the 
roughly 7- Pi cycle. 

The bottom Irf of Fig. [3] was computed from the flrst 
600 pulses of the 610-MHz observation (using overlapping 
256-length ffts). It shows both a low frequency feature and 
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Figure 5. The top plot shows the correlation coefficient of pulses 
800-1100 of the 325-MHz PS in Figs. [3] & fusing only the trail- 
ing longitude regions. The bottom display gives the longitude- 
longitude correlation for a delay of 147 Pi over a longer PS of 
800-1700. The main panel gives the correlation according to the 
color bar on the right, and the average profiles are plotted on 
both axes for reference. The ordinate has been delayed with re- 
spect to the abscissa. The 7-Pi cycle is clearly seen in top plot. 
Note, however, the extraordinary peak at 147 Pi which shows 
some 30% correlation in the top plot and is clearly seen in the 
trailing component in the bottom display. 



two 'sidelobes' adjacent to the primary modulation feature, 
possibly indicating a tertiary modulation as in B0943+10 
(see DROl). The primary modulation has a period of some 
7.11±0.01 Pi, whereas the entire PS gives 7.34±0.01 Pi, sug- 
gesting that the primary modulation frequency is somewhat 
variable. The low frequency feature can be measured over 
an interval of up to 1000 pulses, and the hrf (not shown) in- 
dicates that it represents a mixture of amplitude and phase 
modulation with a period of about 145 Pi. The two 'side- 
bands' are not equally spaced, but the average interval is 
about 0.013 c/Pi. 



VII. EVIDENCE OF A TERTIARY 
MODULATION CYCLE 

We now enquire whether the 7- Pi cycle is a part of a longer 
tertiary modulation cycle as might be expected if the ob- 
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Figure 6. Correlation plots as in Fig. [5] of pulses 800-1900 of 
the 610-MHz PS using both the leading and trailing components. 
Note the peak at 147 Pi . High correlation regions associated with 
both of these components are clearly visible in the bottom display 
at a delay of 147 Pi (see text for details). 



served modulation were produced by a rotating-subbeam 
"carousel" as found, in particular, for B0943+10 (DR99, 
DROl). The long period Irf feature in the 610-MHz PS 
(Fig. O bottom) is suggestive of such a cycle — as are its pair 
of sidebands. It is interesting to note that such a feature 
was also found in recent Westerbork observations at both 
92 (WSE) and 21 cms (WSE)— but its period there can- 
not be determined accurately. However, an autocorrelation- 
function (ACF) analysis of several sections of the the 325- 
MHz, pulse 800-1944 interval identifies a cycle of similar 
length. The top panel of Figure [5] gives the result of this 
analysis done for the outer trailing conal component of the 
800-1100 Pi, where the 7- Pi corrugation is prominent as ex- 
pected. Note the prominent peak at a lag of 147 Pi, which 
is close to 20 times P3. Indeed, this added positive correla- 
tion peak is superposed on the 20th 7- Pi corrugation. This 
peak has an amplitude of about 30%, which is 8 times the 
expected error in the ACF. It is interesting to note the struc- 
ture in the 7- Pi corrugation, which also has a periodicity of 
about 147 Pi. The ACF seems to be consistent with a func- 
tional form of the type cos(27r(i/147) cos(27rd/7), where d is 
the delay in units of Pi. Such orderly variations in the corre- 
lation could reflect non- uniform patterns in either the spac- 
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Figure 7. The top plot shows the cross-correlation between com- 
ponents I and delayed V of pulses 800-1944 of the 325-MHz PS. 
Note the correlation peaks around +94-95 Pi of up to 12%, sig- 
nificantly higher than in the adjacent regions; however, no ex- 
traordinary correlation is seen at a lag of -53 Pi. The +94- Pi 
delay longitude-longitude correlation is shown in the lower dis- 
play, where the correlation between components I and delayed V 
is evident. 



ing or amplitude of the carousal beamlets. The longitude- 
longitude correlation at a lag of 147 Pi (for a slightly longer 
PS of 800-1700) is displayed in the lower part of Fig. [5] This 
sequence exhibits slightly lower (yet significant) correlation 
in the strong trailing conal component. The ACF analysis of 
the 610-MHz observation also shows significant correlation 
at 146 Pi with a peak at 147 Pi as illustrated in the top 
panel of Figure (6] Unlike the 325-MHz PS, here significant 
correlation at 147 Pi is seen for both the leading and the 
traiUng components as patches of red along the diagonal in 
the longitude- longitude correlation map (bottom display). 

If this 147- Pi peak is indicative of the tertiary modula- 
tion, which might be produced by a subbeam carousel sys- 
tem, then we might expect that there should be appreciable 
correlation between components I and V at some part of the 
147- Pi cycle, reflecting the circumstance that the 'beamlets' 
rotate around the outer cone. Consideration of the emission 
geometry (see §111 and Fig.[9]below) suggests that this inter- 
val should be a about 1/3 of the cycle. In both the 325-MHz 
(pulses 800-1944) and 610-MHz PSs (over pulses 1-600; see 
Fig. E] bottom), we find significant cross-correlation between 
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Figure 8. Cross-correlation plots as in Fig. [7] between compo- 
nents I and delayed V of pulses 1-600 of the 610-MHz PS. Again 
we see a prominent peak at a lag of +97 Pi, but little at -50 
Pi. The bottom panel shows the longitude-longitude correlation 
map of the same 610-MHz PS interval in Fig. [3] (bottom) at a 
delay of +97 Pi. Here the PS has been smoothed by 5 samples. 
Note the strong correlation between component I and the delayed 
component V. 



comp. I and delayed comp. V at lags of 94-95 and 96-97 
Pi, respectively, as shown in Figs. [7] and (8] We find no ex- 
traordinary cross-correlation between comp. V and delayed 
comp. I, which we had expected at lags of 53-54 Pi at 325 
MHz (see Fig. E)) and 50-51 Pi at 610 MHz (if we assume 
the exact circulation time is 147-Pi). This might result if 
the actual peak lag fall between two integers. These vari- 
ous correlations reiterate that the carousel 'beamlets' rotate 
positively with longitude. 

Finally, the bugaboo question is whether these re- 
sponses can be aliased, and the answer surely is yes, but it 
is also very unlikely that they are aliases of higher frequency 
responses. Were the primary 7- Pi modulation a first-order 
alias, its true frequency would be 0.86 c/Pi and period cor- 
respondingly 1.16 Pi. Such a modulation would have to be 
very very stable to correlate in the manner we find in Fig. \5\ 
and the 147-Pi peak would be indicative of 126 beamlets! 
We thus appear safe in proceeding on the assumption that 
the 7- and 147- Pi features are not aliased. 
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VIII. CAROUSEL CONFIGURATION 

We are now in a position to compute polar maps correspond- 
ing to the subpulse-modulation patterns discussed above. A 
carousel circulation time P3 of around 147 Pi was indicated 
both directly by the low frequency feature and by ACFs of 
the PSs (Fig. E}. This, together with the 7- Pi P3 feature, 
strongly suggests a carousel with 20 beamlets. We argued 
above that it is highly unlikely that the features are aliased, 
and longitude- longitude correlations appear to confirm that 
the carousel rotates positively with longitude through the 
outer conal components. The longitude of the magnetic axis 
is taken at the midpoint between the outer conal component 
pair, the putative longitude origin used in most of the fig- 
ures above. And while we have not been able to determine a 
and f3 more accurately, the "working" values of 25 and +1.8° 
are sufficiently accurate for our present purposes — the latter 
implying an "outside" or poleward traverse of the sight line. 

Figure [9] displays carousel maps for two sections of the 
325-MHz PS. Both use a local value of P3, computed as 20 
times the P3 determined within the same interval (7.42 and 
7.21 Pi, respectively). The upper carousel map shows the 
average pattern of three rotations during a relatively stable 
interval (pulses 1101-1545; see Fig. [S] top left); whereas the 
bottom map depicts the average configuration during a less 
ordered four-Ps segment (#1-588; see Fig. [3l right). There 
is a good deal to see in these diagrams. First, there is sub- 
stantial irregularity even in the upper map, such that 10 
regularly spaced beamlets are only seen over half the map. 
In the lower map, some beamlets have an azimuthal spacing 
near (360° /20=) 18°, but many do not, so that the broad 
primary modulation feature in the latter figure is not sur- 
prising. Second, the diagrams show virtually no correspon- 
dence between the beamlet patterns in the inner and outer 
conal rings, as expected from our earlier Irf analyses. 



IX. SUMMARY AND CONCLUSIONS 

Five-component, core/double-cone pulsars are relatively 
rare in the normal pulsar population {e.g., Rankin 1993). 
B1237+25 (e.g., Srostlik & Rankin 2005) is by far the best 
studied example, and B 1857-26 has frequently been com- 
pared to it almost as a twin. For the latter, however, little 
beyond average-profile studies have been available until very 
recently, and the present study provides the first sensitive 
PS analysis. 

Like B1237+25, B1857-26 exhibits both strong, regu- 
lar subpulse modulation and null pulses. In the latter, the 
P3 is longer, about 7.4 Pi and is confined to the outer 
conal components. Our two observations have so far identi- 
fied only one behaviour, or profile mode, as against several 
in B1237+25. The nuh fraction is also larger in B1857-26, 
about 20%±3%, and twice the value of 10%=b2.5% reported 
by Ritchings (1976). And also like the former, the nulls in 
this pulsar are short, typically one pulse and no more than 4 
pulses in our observations. Possibly, these are in fact pseudo 
nulls [as in e.g., B0834+06 (Rankin & Wright 2007a) and 
B1133+16 (Herfindal & Rankin 2007)] because the maxi- 
mum null length is about P3/2 — and thus they represent 
"empty" sightline traverses through the subbeam carousel 
system. If so, the 4- Pi nulls may occur preferentially in or- 
derly or disorderly intervals of the observations. 
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Figure 9. Polar maps constructed using the the respective pulse 
1101-1545 (top) and 1-588 (bottom) PSs. Note the roughly reg- 
ular pattern of 20 'beamlets' in the top map and the irregular 
pattern in the lower one. The magnetic axis is at the centre of 
the diagram and the "closer" rotational axis is upward, such that 
the sightline track sweeps through the pattern as indicated. Here 
the star rotates counterclockwise, causing the sightline to cut the 
clockwise-rotating subbeam pattern from right to left; see DROl 
for further details. The side panels give the "base" function (which 
has not been subtracted from the maps here), and the lower pan- 
els show the radial form of the average beam pattern. 



Analyses of B1857-26 (unlike B1237+25) provide evi- 
dence of a tertiary modulation period P3 of some 147 Pi 
at both 325 and 610 MHz. Indications of this tertiary pe- 
riod were first identified via an ACF analysis {e.g., Fig. [S} 
and then corroborated by a low frequency Irf feature (Fig. O 
bottom). While the P3 appears to vary in a range around 
146 to 148 Pi (just as measured values of P3 vary between 
about 7.3-7.4 Pi), we have found that in any given interval 
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the circulation time remains 20 times the primary modu- 
lation period. Numerically, then, some 20 beamlets would 
then nominally comprise the carousel pattern, and the more 
or less regular (or irregular) intervals in our PSs appear to 
reflect differences in the constancy of spacing and possibly 
number of beamlets at a given time. 

The 147 Pi circulation time is again very much larger 
than the value predicted by the R&S (1975) theory. For 
B1857-26, with a 612-ms Pi and a B of 3x10^^ G, this the- 
ory predicts only 5.6 Pi. If ExB forces drive the circulation, 
then it may be that they operate over a larger height range 
than envisioned in the above theory, as suggested, e.g., by 
Hibschman & Arons (2001a,b) and Harding et al (2002). Al- 
ternatively, a vaccum gap partially screened by thermionic 
ion flow in the inner acceleration region above the polar cap 
has been considered (Gil et al 2003) to explain the observed 
slow drift rates in pulsars. 

The PA traverse of B 1857-26 presents an interesting 
case for further analysis. If it is assumed that most of the 
star's emission stems from a single 0PM, then the PA tra- 
verse exceeds 180° by an significant angle. Alternatively, a 
different mode may predominate on the leading and trailing 
sides of the profile, or conceivably, the core emission may 
not follow the RVM. The overall symmetry of the profile, 
however, does not necessarily support a change in the dom- 
inant 0PM, as we know of no other good example of a pul- 
sar with such a marked 0PM asymmetry about the profile 
centre. Such asymmetry is sometimes seen in profiles that 
represent an oblique sight line traverse, but we know of no 
well vetted example of such asymmetry where the sightline 
cuts the beam pattern centrally. 

Finally, the central core component of PSR B 1857-26 
provided us little basis for further analysis. Its shape is not 
at all Gaussian, but otherwise it appears rather steady in 
both intensity and polarization, showing none of the effects 
exhibited by B0329+54 (Mitra et al 2007). 
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